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Abstract 

A comprehensive theoretical description is given for the spin relaxation 
and diffusion in solids. The formulation is made in a general statistical- 
mechanical way. The method of the nonequilibrium statistical operator 
(NSO) developed by D. N. Zubarev is employed to analyze a relaxation 
dynamics of a spin subsystem. Perturbation of this subsystem in solids 
may produce a nonequilibrium state which is then relaxed to an equilib- 
rium state due to the interaction between the particles or with a thermal 
bath (lattice). The generalized kinetic equations were derived previously 
for a system weakly coupled to a thermal bath to elucidate the nature 
of transport and relaxation processes. In this paper, these results are 
used to describe the relaxation and diffusion of nuclear spins in solids. 
The aim is to formulate a successive and coherent microscopic description 
of the nuclear magnetic relaxation and diffusion in solids. The nuclear 
spin-lattice relaxation is considered and the Gorter relation is derived. 
As an example, a theory of spin diffusion of the nuclear magnetic mo- 
ment in dilute alloys ( like Cu-Mn) is developed. It is shown that due 
to the dipolar interaction between host nuclear spins and impurity spins, 
a nonuniform distribution in the host nuclear spin system will occur and 
consequently the macroscopic relaxation time will be strongly determined 
by the spin diffusion. The explicit expressions for the relaxation time in 
certain physically relevant cases are given. 

PACS numbers: 05.30.-d, 05.60.-k, 76.20.+q, 76.60.-k 



1 INTRODUCTION 

For many years there has been considerable interest, experimental and theo- 
retical, in relaxation processes occurring in various spin systems, especially the 
nuclear spin systems in solids and liquids [U EI El El El El El El El HDl HH H21 H31 HH 
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1371138) . In ordinary spin resonance experiments, spins are subject to an applied 
magnetic field ho and make a precessional motion around it. Local fields pro- 
duced by interactions of the spins with their environments act as relatively weak 
perturbations to the unperturbed precessional motion. In quantum-mechanical 
language, the external field gives rise to the Zeeman levels for each spin and the 
interactions are perturbations to these quantum states. In a nuclear-magnetic 
resonance (NMR) experiment, the nuclear spin system absorbs energy from the 
externally applied radio-frequency field and transfers it to the thermal bath 
or reservoir provided by the lattice through the spin lattice interaction. The 
coupled nuclear spins in a solid with very slow spin-lattice relaxation time T\ 
comprise a quasi-isolated system which for many purposes can be treated by 
thermodynamic methods. The spin-spin relaxation time is denoted by T 2 . The 
other system, called the lattice, contains all other degrees of freedom, phonons, 
translational motion of conduction electrons, etc. It is at a temperature T that 
it is considered stable. A macroscopic approach to the description of magnetic 
relaxation was proposed by Bloch 3! . He proposed a phenomenological equation 
describing the motion of nuclear-spin system subjected to both a static and a 
time varying magnetic field 

dM - M x -» M y ~ M -M zr 

— =1 M xh- —i - —±3 + k 

at T 2 T 2 Ti 

where the external field h is taken to be of the form h = h$k + 2hi(t)costuti. 
This equation successfully describes a wide variety of magnetic resonance exper- 
iments, although to obtain a valid description of low-frequency phenomena, it is 
necessary to modify the original equation so that relaxation takes place toward 
the instantaneous magnetic field. In an NMR experiment, the absorption of 
energy from the applied rf field produces either an increase in the energy of the 
spin system or a transfer of energy from the spin system to the lattice. The 
latter process requires a time interval of the order of spin-lattice relaxation time 
T\. The characteristic time T 2 determines the relaxation of the transversal spin 
components due to the spin-spin interactions. 

The relaxation processes in spin systems have been investigated by a number 
of authors |Hlinil2nillIll221ISll23l2Sll23l2HlE3ISZIto obtain qualitative and 
quantitative information about irreversible spin-spin and spin-lattice processes 
in spin systems. The method of many of these papers was to develop an equa- 
tion of motion for the reduced density matrix |12l 1131 [TBI 117) describing the spin 
system, and was found to be most useful when the perturbation responsible 
for the relaxation of the spin system had a very short correlation time. In the 
equation-of-motion approach, the specification of the initial conditions involves 
the assumption of some explicit form for the density matrix describing the sys- 
tem (the system includes both the spin and its surroundings, which in the case 
studied below will be the conduction electrons in a metal). This problem is 
very attractive from the point of view of irreversible statistical mechanics since 
a general model of magnetic resonance consists of a driven system of interest 
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in interaction with a heat bath. Stcnholm and ter Haar JS] have analyzed the 
basic assumptions which are necessary for statistical-mechanical derivation of 
the Bloch equation and the role of the thermal bath. 

An important concept in the interpretation of spin-lattice relaxation phenomena 
was provided by the thermodynamic theory of Casimir and Du Pre pQ. They 
considered the magnetic crystal to be composed of two subsystems which could 
be assigned two different temperatures. One subsystem contained the magnetic 
degrees of freedom. The other subsystem, called the lattice, contains all other 
degrees of freedom. Then, the idea of spin temperature was extended and sev- 
eral distinct temperatures for magnetic subsystems (Zeeman, dipole-dipole, etc. 
) were introduced 122]. (Note, however, some special exclusions In gen- 

eral, the state of the total system to be composed of a few subsystems may be 
described approximately by a density matrix of the form 

p ~ exp[-(#i/feTi) - (H 2 /kT 2 ) - (H 3 /kT 3 ) . . .], 

with a number of quasi-invariant energies TV(iij) and a number of distribution 
parameters T i ~ 1 . Nuclear relaxation in weak applied fields was first treated by 
Redfield |5] and by Hebel and Slichter JU|, using the idea of spin temperature. 
Redfield theory is the semiclassical density operator theory of spin relaxation. 
It was Bloembergen [5J who first formulated that the magnetization of spins in 
a rigid lattice could be spatially transported by means of the mutual flipping of 
neighboring spins due to dipole-dipole interaction. This idea permitted one to 
explain the significant influence of a small concentration of paramagnetic impu- 
rities on spin-lattice relaxation in ionic crystals. He used a quantum-mechanical 
treatment ( first-order perturbation theory) and showed that the transport equa- 
tion for magnetization was a diffusion equation. In this simple approximation 
he calculated the diffusion constant D. In other words, we can roughly represent 
the relaxation dynamics as 



= -A{W<f)) - TO>°] + £>(f)V 2 (/*(r)) 
d(I z (r)) (I z (r)) - (r(r)}° 



(1) 



dt Ti 



1 I 1 

% K lf I + Tf 

where I z is the z-component of the nuclear spin operator. 

Since then, many authors have formulated the general theory of the spin re- 
laxation processes in solids from the standpoint of statistical mechanics or ir- 
reversible thermodynamics. An improvement in the general formulation of the 
theory was achieved by Kubo and Tomita [S] in their treatment of magnetic 
resonance absorption via a linear theory of irreversible processes. In this the- 
ory the important quantities are frequency-dependent susceptibilities which are 
expressed in terms of spin correlation functions. Buishvili |40| developed a 
quantum-statistical theory of the dynamic polarization of nuclei by taking into 
account diffusion of nuclear spins as well as dipole interaction of electron spins. 
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Buishvili and Zubarev (23] developed a successive theory of spin diffusion in 
crystals. The nuclear diffusion in diamagnetic solids with paramagnetic impu- 
rities was analyzed by the method of the statistical operator for nonequilib- 
rium systems. The Bloembergen equation jS], whose coefficients are explicitly 
expressed through certain correlation functions, was obtained. The theory of 
nuclear spin diffusion in the ferromagnets of certain type was considered in 
Ref. 0J. The theory of the dynamic polarization of nuclei and nuclear relax- 
ation for the case of strong saturation was analyzed in Ref. 02] • The influence 
of a strong NMR saturation on spin diffusion was considered in Ref. 03| . The 
time of spin-lattice relaxation was calculated for nuclei when spin diffusion was 
taken into account under conditions of a strong NMR saturation. The influence 
of exchange interactions between nuclear spins on the dynamic polarization of 
nuclei was considered in Ref. 01]. Borckmans and Walgraef 05] formulated a 
theory of Zceman and dipolar energy diffusion in paramagnetic spin systems 
in the frame the general theory of irreversible processes developed by Prigogine 
and co-workers. Buishvili and Giorgadze 46 investigated general theory of spin 
diffusion within the nonequilibrium statistical operator approach. A consistent 
quantum statistical investigation of saturation of a nonuniformly broadened 
EPR line was carried out in Ref. 07] by taking into account spectral diffusion 
and the dipole-dipole reservoir. Role of the flip-flip and flip-flop transitions for 
the dynamic polarization of nuclei was analyzed in Ref. |48|. The theory of 
spin-lattice relaxation in crystals with paramagnetic impurities was discussed 
in paper of Bendiashvili, Buishvili and Zviadadze 0^]- The analysis of the 
role of the interaction between a few subsystems for the construction of the 
nonequilibrium density matrix was discussed from a general point of view by 
Buishvili and Zviadadze J2QJ. The application of the nonequilibrium statistical 
operator method to the case of relaxation in dilute alloys has been considered 
by Fazleev j^|S2]. The influence of relation between thermal capacities of nu- 
clear spin subsystem and the reservoir of electron spin-spin interaction on spin 
kinetics, especially in the low-temperature case when the spin polarization of 
subsystems are high enough was analyzed in detail by Tayurskii |5.'i| for the case 
of insulators. 

Robertson |27| derived an equation of motion for the total magnetic moment 
of a system containing a single species of nuclear spins in an arbitrarily time- 
dependent external magnetic field. He derived a generalization of Bloch's phe- 
nomenological equation for a magnetic resonance. In papers j^|22], the general 
quantum-statistical-mechanical approach to the problem of spin resonance and 
relaxation which utilized a projection operator technique was developed. From 
the Liouville equation for the combined system of the spin subsystem and the 
thermal bath a non-Markoffian equation for the time development of the statis- 
tical density operator for the spin system alone was derived. The memory effects 
were taken into account in the application of the method of the statistical oper- 
ator for nonequilibrium systems to magnetic relaxation problem by Nigmatullin 
and Tayurskii |54| . Romero-Ronchin, Orsky and Oppcnheim |55| used a pro- 
jection operator technique for derivation of the Redfield equations |8]. In their 
paper the relaxation properties of a spin system weakly coupled to lattice 
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degrees of freedom were described using an equation of motion for the spin den- 
sity matrix. This equation was derived using a general weak coupling theory 
which was previously developed. To second order in the weak coupling parame- 
ter, the results are in agreement with those obtained by Bloch, Wangsness and 
Redfield, but the derivation does not make use of second order perturbation 
theory for short times. The authors claim that the derivation can be extended 
beyond second order and ensures that the spin density matrix relaxes to its ex- 
act equilibrium form to the appropriate order in the weak coupling parameter. 
In this work, we present a complementary theory which examines the relaxation 
dynamics of a spin system in the approach of the noncquilibrium statistical op- 
erator. It uses a general formalism from a previous study for a system that is in 
contact with a thermal bath ( a "lattice" ) and relax to the equilibrium state. 
The aim of this paper is to show how the general theory of irreversible processes 
allows a theoretical study of such phenomena without postulated equations of 
phenomenological assumptions. One of our purposes in this paper is to present 
a unified statistical mechanical treatment of spin relaxation and spin diffusion 
phenomena. The transport of nuclear spin energy in a lattice of paramagnetic 
spins with magnetic dipolar interaction plays an important role in many re- 
laxation processes. In this paper the microscopic derivation of an expression 
for the longitudinal relaxation time of bulk metal nuclear spins by dilute local 
moments is performed taking into account spin diffusion processes within the 
nonequilibrium statistical operator approach. 

In the next section, we establish the notation and briefly present the main ideas 
of the NSO approach. This section includes a short summary of the derivation 
of the generalized kinetic and rate equations with the NSO method. Section 
2 serves as an extended introduction to the present paper. In section 3, the 
dynamics of the nuclear spin system is analyzed. We consider the application 
of the established equations to the derivation of the relaxation equations for 
spin systems. Special attention is given to the problem of spin relaxation and 
diffusion in section 4. The case of nuclear spin diffusion in dilute magnetic al- 
loys is discussed in some detail in section 4.2. The final section contains some 
concluding remarks concerning the results obtained. 

2 BASIC NOTIONS 

The statistical mechanics of irreversible processes in solids, liquids, and complex 
materials like a soft matter are at the present time of much interest jSSl E7| 1581 
159) . The central problem of noncquilibrium statistical mechanics is to derive a 
set of equations which describe irreversible processes from the reversible equa- 
tions of motion [571159) . The consistent calculation of transport coefficients is of 
particular interest because one can get information on the microscopic structure 
of the condensed matter. During the last decades, a number of schemes have 
been concerned with a more general and consistent approach to transport the- 
ory |57ll60ll5Tll62| . These approaches, each in its own way, lead us to substantial 
advances in the understanding of the nonequilibrium behavior of many-particle 
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classical and quantum systems. This field is very active and there are many as- 
pects to the problem Our purpose here is to discuss the derivation, within 
the formalism of the nonequilibrium statistical operator 60 , of the generalized 
transport and kinetic equations. On this basis we have derived, by statistical 
mechanics methods, the kinetic equations for a system weakly coupled to a ther- 
mal bath (53] . Our motivation for presenting this alternative derivation is based 
on the conviction that the NSO method provides some advantages in displaying 
the physics of the relaxation processes. 

2.1 Outline of the Nonequilibrium Statistical Operator 
Method 

In this section, we briefly recapitulate the main ideas of the nonequilibrium 
statistical operator approach [HO] E3] for the sake of a self-contained formula- 
tion. The precise definition of the nonequilibrium state is quite difficult and 
complicated, and is not uniquely specified. Since it is virtually impossible and 
impractical to try to describe in detail the state of a complex macroscopic sys- 
tem in the nonequilibrium state, the method of reducing the number of relevant 
variables was widely used. A large and important class of transport processes 
can reasonably be modelled in terms of a reduced number of macroscopic rele- 
vant variables. There are different time scales and different sets of the relevant 
variables |651 166| . e.g. hydrodynamic, kinetic, etc. The most satisfactory and 
workable approach to the construction of Gibbs-type ensembles for the nonequi- 
librium systems, as it appears to the writer, is the method of nonequilibrium 
statistical operator (NSO) developed by D. N. Zubarev (HDJ - The NSO method 
permits one to generalize the Gibbs ensemble method to the nonequilibrium 
case and construct a nonequilibrium statistical operator which enables one to 
obtain the transport equations and calculate the kinetic coefficients in terms of 
correlation functions, and which, in the case of equilibrium, goes over to the 
Gibbs distribution. Although this method is well known, we shall briefly recall 
it, mostly in order to introduce the notation needed in the following. 
The NSO method sets out as follows. The irreversible processes which can 
be considered as a reaction of a system on mechanical perturbations can be 
analyzed by means of the method of linear reaction on the external perturba- 
tion [HH]- However, there is also a class of irreversible processes induced by 
thermal perturbations due to the internal inhomogeneity of a system. Among 
them we have, e.g., diffusion, thermal conductivity, and viscosity. In certain 
approximate schemes it is possible to express such processes by mechanical per- 
turbations which artificially induce similar nonequilibrium processes. However, 
the fact is that the division of perturbations into mechanical and thermal ones 
is reasonable in the linear approximation only. In the higher approximations in 
the perturbation, mechanical perturbations can effectively lead to the appear- 
ance of thermal perturbations. 

The NSO method permits one to formulate a workable scheme for description 
of the statistical mechanics of irreversible processes which include the thermal 
perturbation in a unified and coherent fashion. To perform this, it is necessary 
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to construct statistical ensembles representing the macroscopic conditions de- 
termining the system. Such a formulation is quite reasonable if we consider our 
system for a suitable large time. For these large times the particular properties 
of the initial state of the system are irrelevant and the relevant number of vari- 
ables necessary for description of the system reduces substantially |rI5] . 
The basic hypothesis is that after small time-interval r the nonequilibrium distri- 
bution is established. Moreover, it is supposed that it is weakly time-dependent 
by means of its parameter only. Then the statistical operator p for t > t can 
be considered as an "integral of motion" of the quantum Liouvillc equation 

Here ^ denotes time differentiation with respect to the time variable on which 
the relevant parameters F m depend. It is important to note once again that 
p depends on t by means of F m (t) only. We may consider that the system is 
in thermal, material, and mechanical contact with a combination of thermal 
baths and reservoirs maintaining the given distribution of parameters F m . For 
example, it can be the densities of energy, momentum, and particle number 
for the system which is macroscopically defined by given fields of temperature, 
chemical potential and velocity. It is assumed that the chosen set of parameters 
is sufficient to characterize macroscopically the state of the system. The set of 
the relevant parameters are dictated by the external conditions for the system 
under consideration and, therefore, the term appears as the result of the 
external influence upon the system. Due to this influence precisely, the behavior 
of the system is nonstationary. 

In order to describe the nonequilibrium process, it is necessary also to choose 
the reduced set of relevant operators P m , where m is the index ( continuous 
or discrete). In the quantum case, all operators are considered to be in the 
Heisenberg representation 

P m (t) = exp(^)P m exp(^^) (3) 

where H does not depend on the time. The relevant operators may be scalars 
or vectors. The equations of motions for P m will lead to the suitable "evolution 
equations" |60| . In the quantum case 

2M)-l [Pm ( t ,,*]=0. (4) 

The time argument of the operator P m (t) denotes the Heisenberg representation 
with the Hamiltonian H independent of time. Then we suppose that the state 
of the ensemble is described by a nonequilibrium statistical operator which is a 
functional of P m (t) 

p(t)=p{...P m (t)...} (5) 
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Then p(t) satisfies the Liouville equation (J2J. Hence the quasi-equilibrium ( 
" local-equilibrium" ) Gibbs-type distribution will have the form 



Pq = Q- 1 exp (- F m (t)P m J (6) 

where the parameters F m (t) have the meaning of time-dependent thermody- 
namic parameters, e.g., of temperature, chemical potential, and velocity ( for 
the hydrodynamic stage), or the occupation numbers of one-particle states (for 
the kinetic stage). The statistical functional Q q is defined by demanding that 
the operator p q be normalized and equal to 

Q q = Tr exp ^- F m (t)P rr )j (7) 

This description is still very simplified. There are various effects which can 
make the picture more complicated. The quasi-equilibrium distribution is not 
necessarily close to the stationary stable state. There exists another, com- 
pletely independent method for choosing a suitable quasi-equilibrium distribu- 
tion |56l 1571 IrTTl 167) . For the state with the extremal value of the informational 
entropy [371 I^T] 

S = -Tr(p]np), (8) 

provided that 

Tr(pP m ) = (P m ) q ; Trp=l, (9) 

it is possible to construct a suitable quasi-equilibrium ensemble. Then the cor- 
responding quasi-equilibrium ( or local equilibrium ) distribution has the form 

Pq = exp F m (t)P m \ = exp(-S(t, 0)) (10) 

n = In Tr exp ^- ^ F m (t)P m ^J , 

where S(t, 0) can be called the entropy operator. The form of the quasi- 
equilibrium statistical operator was constructed so as to ensure that the ther- 
modynamic equalities for the relevant parameters F m (t) 

SlnQ q _ SQ _ SS 

6F m (t) ~ 5F m (t) ~ {rm)q > 5(P m ) g -* m[t) [n) 

are satisfied. It is clear that the variables F m (t) and (P m )q are thermodynami- 
cally conjugate. Here the notation used is (. . .) q = Tr(p q . . .). 
By definition a special set of operators should be constructed which depends on 
the time through the parameters F m (t) by taking the invariant part of the oper- 
ators F m (t)P m occurring in the logarithm of the quasi-equilibrium distribution, 
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i.e., 



r° 

B m (t) = F m {t)P m = e / e etl F m {t + t^P^dh = (12) 

J — OG 

F m (t)P m - f dhe etl (F m (t + ti)P m (ti) + F m (t + ii)P m (ii)) 



where (e — > 0) and 



1 ' rfi? m(*) 

Pin — TT" I -f?n , -ff I , Pm{t) — T • 

l« at 



The parameter e > will be set equal to zero, but only after the thermodynamic 
limit has been taken. Thus, the invariant part is taken with respect to the 
motion with Hamiltonian H. The operation of taking the invariant part, of 
smoothing the oscillating terms, is used in the formal theory of scattering 
to set the boundary conditions which exclude the advanced solutions of the 
Schrodinger equation. The P m (t) will be called the integrals ( or quasi- integrals 
) of motion, although they are conserved only in the limit (e — ► 0). It is clear that 
for the Schrodinger equation such a procedure excludes the advanced solutions 
by choosing the initial conditions. In the present context this procedure leads 
to the selection of the retarded solutions of the Liouville equation. 
It should be noted that the same calculations can also be made with a deeper 
concept, the methods of quasi-averages 160) . 

- ^[ln Pe ,H] = -e(lnp e \n Pq ), (13) 

where (e — > 0) after the thermodynamic limit. The required nonequilibrium 
statistical operator is defined as 

o 

eti 



Pe =p e (t,Q) = p g (t,0) =e / dhe^p^t + t^h) (14) 

J — oo 

Hence the nonequilibrium statistical operator can then be written in the form 



Q 1 exp ^- ^2 B mj 

Q exp{- F m{t)Pm +J2 f dt ^ £tl [^(* + ^) P m(h) + 

™ ™ J — OG 



p 

A) 



F m {t + tx)P m {tx)]} (15) 

Now we can rewrite the nonequilibrium statistical operator in the following 
useful form: 

o 

eti 



p(t, 0) = exp -e / dt x e etl lnp 9 (i + t\,t x ) 



exp(lnp,(i,0)) = exp(-S(t,0)) (16) 
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The average value of any dynamic variable A is given by 

(A) = lim Tr(p(t,0)A) (17) 

£—►0 + 

and is, in fact, the quasi-average. The normalization of the quasi-equilibrium 
distribution p q will persists after taking the invariant part if the following con- 
ditions are required 

Tr(p(t, 0)P m ) = (P m ) = (P m ) q ; Trp = 1 (18) 

Before closing this section, we shall mention some modification of the "canon- 
ical" NSO method which was proposed in Ref. [S7| and which one has to take 
into account in a more accurate treatment of transport processes. 



2.2 Transport and Kinetic Equations 

It is well known that the kinetic equations are of great interest in the theory of 
transport processes. Indeed, as it was shown in the preceding section, the main 
quantities involved are the following thermodynamically conjugate values: 

The generalized transport equations which describe the time evolution of vari- 
ables (P m ) and F m follow from the equation of motion for the P m , averaged 
with the nonequilibrium statistical operator Ijlfcifl . It reads 

^--Eidssbj*"'* ^ (() = ?»«™ <20) 

The entropy production has the form 

S(t) = (5(t,Q)) =-J2{Pm)F m (t) = -Y, tF MF 7i) P n(t)F m (t) (21) 

m n,m m w n \ / 

These equations are the mutually conjugate and with Eq. I|19fl form a complete 
system of equations for the calculation of values (P m ) and F m . 



2.3 System in Thermal Bath: Generalized Kinetic Equa- 
tions 

In paper |64j we derived the generalized kinetic equations for the system weakly 
coupled to a thermal bath. Examples of such systems can be an atomic ( 
or molecular) system interacting with the electromagnetic held it generates as 
with a thermal bath, a system of nuclear or electronic spins interacting with 
the lattice, etc. The aim was to describe the relaxation processes in two weakly 
interacting subsystems, one of which is in the nonequilibrium state and the 
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other is considered as a thermal bath. The concept of thermal bath or heat 
reservoir, i.e., a system that has effectively an infinite number of degrees of 
freedom, was not formulated precisely. A standard definition of the thermal bath 
is a heat reservoir defining a temperature of the system environment. From a 
mathematical point of view , a heat bath is something that gives a stochastic 
influence on the system under consideration. In this sense, the generalized 
master equation [70] is a tool for extracting the dynamics of a subsystem of 
a larger system by the use of a special projection techniques [7T] or special 
expansion technique |72| . The problem of a small system weakly interacting 
with a heat reservoir has various aspects. Basic to the derivation of a transport 
equation for a small system weakly interacting with a heat bath is a proper 
introduction of model assumptions. We are interested here in the problem of 
derivation of the kinetic equations for a certain set of average values ( occupation 
numbers, spins, etc.) which characterize the noncquilibrium state of the system. 
The Hamiltonian of the total system is taken in the following form: 



Here Hi is the Hamiltonian of the small subsystem, and a\ and a a are the 
creation and annihilation second quantized operators of quasiparticles in the 
small subsystem with energies E a , V is the operator of the interaction between 
the small subsystem and the thermal bath, and H 2 is the Hamiltonian of the 
thermal bath which we do not write explicitly. The quantities <fr a p are the 
operators acting on the thermal bath variables. 

We assume that the state of this system is determined completely by the set 
of averages (P a p) — { a a a p) and the state of the thermal bath by (H 2 ), where 
(. . .) denotes the statistical average with the nonequilibrium statistical operator, 
which will be defined below. 

We take the quasi-equilibrium statistical operator p q in the form 



p,(t)=exp(-5(t,0)), S(t,0)=n(t) + J2 p <*pF<*f3{t)+(3H 2 (24) 



Here F a p(i) are the thermodynamic parameters conjugated with P a p, and /3 is 
the reciprocal temperature of the thermal bath. All the operators are considered 
in the Heisenberg representation. The nonequilibrium statistical operator has 
the form 



H = H 1 +H 2 + V, 



(22) 



where 




(23) 



n = In Tr exp(- £ P a0 F a p{t) - (3H 2 ) 



p(t)=exp(-S(t,0)) 



(25) 
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The parameters F a p(t) are determined from the condition (P a p) = (P a p)q- 
In the derivation of the kinetic equations we use the perturbation theory in a 
"weakness of interaction" and assume that the equality (^ a p)q = holds, while 
other terms can be added to the renormalized energy of the subsystem. For 
further considerations it is convenient to rewrite p q as 

P q = P1P1 = Qq 1 exp(-Lo(*)), (26) 

where 

Pi = Qi 1 exp I - ^2P a pF aP (t) J ; Q 1 =Trexpl-^2P aP F a p(t)\ (27) 

\ a/3 / \ a/3 / 

P2 = Q^e-^; Q 2 = Trexp(-0H 2 ) (28) 
Q q = QiQ2\ L Q =Y,P al 3F a f]{t)+l3H 2 (29) 

a/3 

We now turn to the derivation of the kinetic equations. The starting point is 
the kinetic equations in the following implicit form: 

= U[P aP , H]) = \(E P - E a )(P af} ) + i<[P Q/J , V]) (30) 
at in in in 

We restrict ourselves to the second-order in powers of V in calculating the r.h.s. 
of H30JI ■ Finally we obtain the kinetic equations for (Pap) in the form |M] 



d(P af: 

dt 



2-{E p -E a ){P aP )- i / dtie rfl ([[^.11.n<i)]> 9 (31) 
w n J_oo 



The last term of the right-hand side of Eq. 131fl can be called the generalized " col- 
lision integral" . Thus, we can see that the collision term for the system weakly 
coupled to the thermal bath has a convenient form of the double commutator 
as for the generalized kinetic equations |73) for the system with small interac- 
tion. It should be emphasized that the assumption about the model form of the 
Hamiltonian (12211 is nonessential for the above derivation. We can start again 
with the Hamiltonian (1221) in which we shall not specify the explicit form of Hi 
and V. We assume that the state of the nonequilibrium system is characterized 
completely by some set of average values (Pk) and the state of the thermal bath 
by (H?). We confine ourselves to such systems for which [Hi,Pk] = ^^kiPi- 
Then we assume that (V) q ~ 0, where (. . .) q denotes the statistical average with 
the quasi-equilibrium statistical operator of the form 



Pq 



= Q- 1 exp - £ P k F k (t) - (3H 2 (32) 



and Fk(t) are the parameters conjugated with (Pk)- Following the method 
used above in the derivation of equation (|31|l . we can obtain the generalized 
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kinetic equations for (Pk) with an accuracy up to terms which are quadratic in 
interaction 

d(P k ) i ' '•" 



dt H 



£c fcI (fi)--^ /" dhe^m.VlVih)}), (33) 



Hence l|31[) is fulfilled for the general form of the Hamiltonian of a small system 
weakly coupled to a thermal bath. 

2.4 System in Thermal Bath: Rate and Master Equations 

In section 12731 we have described the kinetic equations for (Pap) in the general 
form. Let us write down equations (|31|l in an explicit form. We rewrite the 
kinetic equations for (Pap) as 

]T {K pv {P av ) + K av {P v p)) + K ap ^(P^) (34) 

V 

The following notation were used 

1 f° 

-]T/ dt 1 e^(^p IM <f >in/ (t 1 )) g = 



1 "° 



2vr 



dtxe^^^pu^q + (^ a (h)^p u )g) 



it) 

1 /" +0 ° . . , . / I 1 



du>Jp v ^ a (uj) 



27r J.oo \ Huj — Ep + E v + is fiuj — E a — E^ — is / 

= K a p^ v (36) 

Let us now remind jSU] that the correlation functions (AB(t)) and (A(t)B) can 
be expressed via their spectral weight function ( or spectral intensity) J(uj) 

1 r +oc 

F AB (t-t l ) = (A(t)B(t')) = — duJcMMt-t')]J A B(uj) (37) 

^ J-oo 



F BA (t' -t) = (B(t')A(t)) = —J dm exp[iu(f - t)]J BA (u>) (38) 

The correlation functions (Qp^^v^ifjq and ((t>vn(ti)&iict)q are connected with 
their spectral intensities in the following way: 

^ /> + oo Jj\ 

( < bu, v 4>~i& (*)}? = 7T~ / dwJ 7 *,u I /(w)exp[-i(o; :L T )t] (39) 

(4>^(t)^is) q = — / du)J lS ,pv{u) exp[i(w + — -)t] (40) 

^7T /_„ ft 
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The above result is similar in structure to the Redheld equation for the spin 
density matrix |S] when the external time-dependent field is absent. Indeed, the 
Redfield equation of motion for the spin density matrix has the form |S] 

VP' 

Here p aa is the a, a' matrix element of the spin density matrix, uj aa > = [E a — 
E a i)h, where E a is energy of the spin state a and R aa ' pp 1 P^^ is the "relaxation 
matrix" . A sophisticated analysis and derivation of the Redfield equation for 
the density of a spin system immersed in a thermal bath was given in |74|. A 
brief discussion of the derivation of the Redfield-type equations in an external 
field is given in Appendix A. 

Returning to Eq. Ij34(l , it is easy to see that if one confines himself to the diagonal 
averages (P aa ) only, this equation may be transformed to give 

d(P aa ) _ (P )-(K )(P ) (41) 



Kaa,Pf3 = ~~2 ^ 0tf3,f3a ( ~~ T ) = Wfj^ a (42) 

h n 

K aa + Kt a = ±J2 Jf3<*M Ef3 Z Ea ) = ^ (43) 

Here Wp^ a and W a ->p are the transition probabilities expressed in the spectral 
intensity terms. Using the properties of the spectral intensities [HO], it is possible 
to verify that the transition probabilities satisfy the relation of the detailed 
balance 

Wp^a _ exp(-/3E a ) 



Finally, we have 

d(Paa 



(44) 



dt 



^W v ^ a {P m ) -J2 w ^AP aa ) (45) 



This equation has the usual form of the Pauli master equation. 

It is well known that "the master equation is an ordinary differential equation 

describing the reduced evolution of the system obtained from the full Heisenberg 

evolution by taking the partial expectation with respect to the vacuum state of 

the reservoirs degrees of freedom" . The rigorous mathematical derivation of the 

generalized master equation |7()l 1711 172*] is rather a complicated mathematical 

problem. 
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3 DYNAMICS OF NUCLEAR SPIN SYSTEM 

In nuclear magnetic resonance one has a system of nuclei with magnetic moment 
p and spins / which are placed in a magnetic field ho- The magnetic moment p 
and momentum of nuclei J ~ hi are related as p = j n J = "f n hl = g n r)l , where 
7„ is the gyromagnctic nuclear factor, g n is the nuclear spectroscopic factor, and 
77 = eh/2Mc is the nuclear magneton. If the spins are otherwise independent, 
their interaction with the imposed field produces a set of degenerate energy 
levels which for a system of N spins are (21 + 1) N in number with the energy 
spacing hui n = /j.ho/1. It should be noted that the method of NMR is most 
powerful and useful in diamagnetic materials. Metals may be studied, although 
there are some technical specific problems. 

In a nuclear-magnetic-resonance experiment, the nuclear spin system absorbs 
energy from the externally applied radio-frequency field and transfers it to the 
thermal bath or reservoir provided by the lattice through the spin-lattice inter- 
actions. The latter process requires a time interval of the order of the spin-lattice 
relaxation time T\. The term "lattice" is used here to denote the equilibrium 
heat reservoir with temperature T associated with all degrees of freedom of the 
system other than those associated with the nuclear spins. 
A great advantage of magnetic resonance method is that the nuclear spin sys- 
tem is only very weakly coupled to the other degrees of freedom of the complex 
system in which it resides and its thermal capacity is extremely small. It is, 
therefore, possible to cause the nuclear spin system itself to depart severely from 
thermal equilibrium while leaving the rest of the material essentially in thermal 
equilibrium. As a consequence, the disturbance of the system other than the 
nuclear spins could be ignored. 

If the nuclei are in thermodynamic equilibrium with the material at temperature 
T in a field ho, a nuclear paramagnetic moment M is produced in the direction 
of ho given by the Curie formula M /h — n^i 2 /3kT, n is the number of nuclei 
per unit volume. 

We can evidently disturb the system from equilibrium by applying radiation 
from outside with quanta of size hw n and with suitable polarization. If the 
equilibrium distribution is disturbed and the population changed, the magneti- 
zation in the z-direction, M z , is different from Mo, say M z . If then left alone, 
M z reverts to M and usually does so exponentially with time, i.e. 

M x {t) = M - (Mo - M h z ) exp{-^} 

The last expression serves to define the spin-lattice relaxation time, T\, and is 
so called because the process involves exchange of magnetic orientation energy 
with thermal energy of other degrees of freedom ( known conventionally as a 
lattice ). All the interactions with the nucleus may contribute to the relaxation 
process so we must add all contributions to 1/T\ 

1 1 11 

J-i J- la J- lj 
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where various contributions to relaxation due to various interactions have been 
added. The relaxation rates may be dominated by one or more different physical 
interactions, so that the observable power spectrum may be the Fourier trans- 
form of functions involving dipole-dipole correlations, electric field gradient- 
nuclear quadrupole moment correlations, etc. 

The dipole-dipole interaction Hamiltonian H ( id between the magnetic moments 
of nuclei may contribute significantly to the nuclear magnetic relaxation pro- 
cess |75) . Consider an explicit interaction between the moments fix and fi2 which 
are distant by r*i2 from each other. Then the interaction is written as 

3(/Zifi 2 )(M2fi2) !A( .\ 
H dd = -3 -g = (46) 

'12 '12 

[Att 1 

~\ — — [2/xf/i|l2,0 - (Ml>2 + [hV-t) Y 2.0 + 

V a t-12 

v/3(mJ M2 + Mi/4)*2,-i + V3(^n z 2 + M1M2 )^2,i 

+ V&nfntY2,-2 + V^fJ- 1 ^J-2 Y 2,2] 

where /i ± = (/i x ±/i y )/v2 and Y 2 , m denote the normalized spherical harmonics 
of the second degree expressed in the form 



15 / 15 

*2,±2 = \l sin 2 6»i2 exp(±2i(/>); Y 2±1 = \ 3- sin 6*12 cos0i 2 exp(±i<j 



r2 ' o = Vifc (3cos2 612 _1) 

The dipole-dipole coupling provides the dissipation mechanisms in the spin sys- 
tem. It acts as time dependent perturbations on the Zeeman energy levels, 
which results in the relaxation of the nuclear magnetization. 
Thus, such a spin system can be described as a superposition of a number of 
subsystems. They are the Zeeman subsystem for each spin species and the 
dipole-dipole subsystem. A weak applied rf field can be considered as an ad- 
ditional subsystem. The coupling inside each subsystem is strong, whereas the 
coupling between subsystems is weak. As a consequence, the subsystems reach 
internal thermal equilibrium independently of each other and one can ascribe a 
temperature, an energy, an entropy, etc., to each of them. Let us note that the 
usual prediction of statistical mechanics that the temperatures of interacting 
subsystems become equal in equilibrium is a direct consequence of the conjec- 
ture that the total energy is the only analytic constant of the motion. 



3.1 The Hierarchy of Time Scales 

A case of considerable practical interest in connection with the phenomenon of 
resonance and relaxation is that of the hierarchy of time scales. In the standard 
situations the interaction between nuclear spins is weak as well as the interac- 
tion with the lattice is weak. As a result, in the NMR case the thermal bath 
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variables change on the fast time scale characterized by ti, c while the spin vari- 
ables change on the slow time scale characterized by T sr . 

First of all, consider the most important concept of spin temperature ■ Actu- 
ally, spin systems are never completely isolated and the concept of spin temper- 
ature is meaningful only if the rate Tq 1 of achievement of internal equilibrium 
is much faster than the spin-lattice relaxation rate T± . For time intermedi- 
ate between To and T sr , the spin temperature exists and can be different from 
the lattice temperature T. The necessary condition for the applicability of spin 
temperature concept is then inequality To ~ T-i <C T±. 

Characteristic times are long in comparison with the time of achievement of in- 
ternal equilibrium in the lattice the but short compared to spin relaxation times 
the < t < T sr . In this case, the second-order perturbation theory is valid in 
the weak spin-lattice coupling parameter. Usually, it is assumed that the time 
the is very short and tq > the- The restriction of ordinary perturbation theory 
generally applied is that it is valid when within the time interval considered the 
density matrix cannot change substantially. Argyres and Kelley |22| removed 
the restriction t^ c < r sr and derived an equation of motion for the spin density 
that depends on the history of the system ,74 ■ 

One of the essential virtues of the NSO method is that it focuses attention, at 
the outset, on the existence of different time scales. Suppose that the Hamilto- 
nian of the spin system can be divided as H = H + V, where H is the dominant 
part, and V is a weak perturbation. The separation of the Hamiltonian into 
Ho and V is not unique and depends on the physical properties of the system 
under consideration. The choice of the operator Ho determines a short time 
scale To- This choice is such that for times t 3> to the nonequilibrium state of 
the system can be described with a reasonable accuracy by the average values 
of some finite set of the operators P m ( I3J) • 

After the short time To, it is supposed that the system can achieve the state 
of an incomplete or quasi-equilibrium state. The main assumption about the 
quasi-equilibrium state is that it is determined completely by the quasi-integrals 
of motion which are the internal parameters of the system. The characteristic 
relaxation time of these internal parameters is much longer than To. Clearly 
then, that even if these quasi-integrals at the initial moment had no defini- 
tive equilibrium values, after the time To, at the quasi-equilibrium state, those 
parameters which altered quickly became the functions of the external param- 
eters and of the quasi-integrals of motion. It is essential that this functional 
connection does not depend on the initial values of the parameters. In other 
words, the operators P m are chosen so that they should satisfy the condition : 
[Phi Hq] — cm.Pi- It is necessary to write down the transport equations ( I20fl 
for this set of "relevant" operators only. The equations of motion for the average 
of other "irrelevant" operators ( other physical variables) will be in some sense 
consequences of these transport equations. As for the "irrelevant" operators 
which do not belong to the reduced set of the "relevant" operators P m , relation 
[Pie, Ho] = J^i c kiPi leads to the infinite chain of operator equalities. For times 
t < To the nonequilibrium averages of these operators oscillate fast, while for 
times t > to they become functions of the average values of the operators. 
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3.2 Nuclear Spin-Lattice Relaxation 

At the earlier stage, the theory of spin relaxation was developed by means of 
quantum mechanical perturbation methods. Here the spin relaxation is studied 
by making use of the method of noncquilibrium statistical operator. We discuss 
in this section an arbitrary nuclear spin system on a lattice in interaction with 
external fields and another system |7S], to be taken eventually to act as a heat 
bath. The bath is considered as a quantum-mechanical system that remained in 
thermodynamic equilibrium while its exchange of energy with the spin system is 
taken into account. We consider the processes occurring after switching off the 
external magnetic field in a nuclear spin subsystem of a crystal. Let us consider 
the behavior of a spin system with the Hamiltonian H n weakly coupled by a 
time-independent perturbation V to a thermal bath (temperature reservoir) or 
a crystal lattice with the Hamiltonian H^. 
The total Hamiltonian has the form 

H = H n +H L + V, (47) 

where 

H n = -a^If; a = lnh (48) 

i 

Here If is the operator of the z-component of the spin at the site i , ho is 
the time-independent external field applied in the z-direction, and 7„ is the 
gyromagnetic coefficient. 

Now we introduce b iX and bi\ the creation and annihilation operators of the spin 
in the site i with the z-component of the spin equal to A , where —I < A < I. 
Then we have 

A A 

and, consequently, 

H n = J2 E Wiy> E x = -a\ (50) 

i\ 

Following section |2~31 we write the Hamiltonian of the interaction as 

v = E E $ ^AtA^ = (5i) 

i it,u 

Here ^iv,i^ are the operators acting only on the "lattice" variables. The term 
" lattice" is used here to denote the equilibrium heat reservoir with temperature 
T associated with all degrees of freedom of the system other than those associ- 
ated with the nuclear spins. Then, in agreement with Ea. pfijl . we construct the 
quasi-equilibrium statistical operator 

P q =PL®Pn, (52) 
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where 

PL = Ql'e-^; Q L = TreM-0H L ) (53) 

„ , sinh ^#a(27 + 1) 

Pn = Q- N exp(-f3 n (t)H n ); Q n = 2 ~ T \ '- (54) 

sinh ^ ' a 

Here f3 n is the reciprocal spin temperature and ./V is the total number of spins 
in the system. 

We now turn to writing down the kinetic equations for average values (rii\) — 
(b iX bi\). We use the kinetic equation in the form 



d(rii 



dt 

where 



J2w^(ii){n iu ) -J2 w x^u(n){n iX }, (55) 



W v ^x{ii) = ^2^, to ^,» A ( — h — ) (56) 
It can be shown that 

(riix) = (n\) = Q" 1 exp[-f3 n E x } 



Wx^, = ^^Wx^u(U); W^x = -^Y,W^x(U) (58) 



Then we obtain 

d(n x ) 
dt 

where 

UX = jv^"— 'w. N . 

i i 

It is easily seen that 

W v ^x = exjp\p{E v - E x )]Wx^ 

Hence, for (3 n we find the equation 

dpn 1 Euxi^ ~ v)Wx^ v (1 ~ exp[-(/3 - !3 n ){Ex - E v )\) exp[-f3 n E x ] 
dt 2 ^ 82 '"Q" 



(59) 



In the derivation of Ea. (|59|l we took into account that (I z ) — v { n u) and 



d{I z ) = l d/? K d 2 lnQ K = lrf/? n 
di a <9/3 2 a 



(((r)2) _ (/ , )2) _ (6Q) 
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In the high-temperature approximation (huj n <C kT )we obtain 

d/3 n _ (3 - P n 
dt Ti ' 

where T\ is the longitudinal time of the spin-lattice relaxation 

The above expression is the well-known Gorter relation 1141 1331 13~5] . 



(61) 



(62) 



4 SPIN DIFFUSION OF NUCLEAR MAGNETIC 
MOMENT 

The concept of spin diffusion was invoked by Bloembergen to explain the 
magnetic relaxation of nuclei in diamagnetic solids, which is due to the in- 
teraction of the nuclear spins with spin of a paramagnetic impurity ion. This 
theoretical approach was further developed in many works |U1 1231 1331 1761 ]77\ . 
In the previous section, we have discussed a simple calculation of the longitudi- 
nal nuclear spin relaxation time within the NSO approach. Here we shall extend 
this treatment in order to obtain a more sophisticated description of the spin 
dynamics. Let us, therefore, work out a general formula, using these ideas. 
Consider a subsystem of interacting nuclear spins I of a crystal which interact 
with the external magnetic field fig and with other subsystems of a crystal. Our 
aim is to derive the evolution equation for the reciprocal spin temperature of 
the Zeeman spin subsystem f3 n (r , t) which is relaxed to the equilibrium after 
switching off the external rf field. The total Hamiltonian has the form 

H = H n + H dd + H L + V, (63) 

where the Zeeman operator H n is given by 

H n = -a ^2 If; a = j n h (64) 

i 

It is convenient to rewrite H n in the following form: 

H n {7) = J2lfh(u n + n i )6(r-n) (65) 

i 

Here f2^ <C 0J n is effective renormalization of the " bare" nuclear spin energy Tko n 
due to the surrounding medium and will be written explicitly below; Hdd is the 
operator of dipole-dipole interaction l|46|l 

H dd = J^M 3(£f)(£f)}, 

ij 
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where r is the distance between the two spins and f = is the unit vector 

in the direction joining them. It was shown [321 155] that the so-called secular 
part of this operator was essential, and in the rest of the paper we will use the 
notation Hdd for the secular part of the operator of dipole-dipole interaction. It 
has the form [521153] 

Here ^ 

Aij = ^(l-3cos 2 %) 

and 9ij is the angle between ho and r7j ; 

Hl is the Hamiltonian of the a thermal bath and V is the operator of interaction 
between the nuclear spins and the lattice. Since our aim is to derive the equation 
for the relaxation of the Zeeman energy, we take the operators H n (f) and Hdd 
as the relevant variables which describe the nonequilibrium state. According to 
the NSO formalism, we now write the entropy operator l|24|) in the form 

S(t, 0) = Q(t) + 0H L + [3 d H d d + J Pn(r, <)W„(r)d 3 r, (67) 

p,(t) = exp(-S(f,0)) 

where (3d and (3 are the reciprocal temperature of dipole-dipole subsystem and 
the thermal bath, respectively. Then, within the formalism of NSO, as de- 
scribed above in section T2. 31 it is possible to derive the corresponding transport 
equations for the nonequilibrium averages (TC n (r)) and (Hdd)- Here we confine 
ourselves to the equation for the (Tt n (r)) since the equations for (3 n (r, t) and (3d 
are decoupled when the external rf field is equal to zero. We need the relations 

= 4 [H n (r) , V] + 4 [H n (r) , H dd ] = K n (r) - divJ{r) (68) 
dt in in 

Here K n (f) is the source term and J(r) is the effective nuclear spin energy 
current 

J{r) = 4 E A U f M(u n + SliW- ri)I+Ii (69) 
Since fij -C u> n , the approximate form of the current is 

f W M ^ E 4w*W" rl)I+I- (70) 

The law of conservation of energy in the differential form can be written as 

(c.f. |231) 

= ~div(J(r)) + (K n (r)) (71) 
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Following the method of calculation of Buishvili and Zubarev we get 



^r^ = - E ^^<f)-^ft*m+ifi»m-fii*w (72) 

fj,i/= 1,2,3 " 

According to Ea.(lfi5|l. we have treated (I z (t)) as a continuum function of spatial 
variables so that when evaluated at the lattice site j, it is equal to (I z (t)}. 
Carrying out a Taylor series expansion |3U] of (I z {t)} about the kth lattice site 
and then evaluating the results at position j yield 

a;®} « vm) + E ^a*(*))ifc^i- + (73) 

a— 1 



5 E 



a 2 



2 V &E a &E0 

a,/3=l 



(I Z (t))\kXkjocX kj l3 + 



where ccfcjc is the a coordinate (a = 1,2,3) in an arbitrary Cartesian coordinate 
system for rfej , and d/dx a (I z (t))\k is the partial derivative of (I z (t)) with respect 
to x a , evaluated at the lattice site k. 

The generalized kinetic coefficients L^ v {f) and L\{f) have the form 



e \ 







L^(r)= / dii / / d 3 g(J A1 (r)exp(-A5(i,0))J,(q,i 1 )exp(A5(i,0))) 9 



J}(r) = f dh f dX f 

J-oc JO J 



(74) 

d 3 q(K n (f) exp(-AS(t, 0))tf n (g, i x ) exp(A5(t, 0))), 

(75) 



The condition (H n {r)) = CH n (r)) q determines the connection of f3 n (r,t) and 
(H n (r)). Equation ( 172) 1 is the diffusion type equation JJ^HZSUHDI- This equation 
describes more fully the local changes of the Zeeman energy due to the relaxation 
and transport processes in the system with the Hamiltonian (|fi3|) . In its general 
form equation ( 172")) is very complicated OEnHEH] and to get a solution, various 
approximate schemes should be used. 

4.1 Evaluation of Spin Diffusion Coefficient 

Let us consider the calculation of the diffusion coefficient. The most obvious 
approximation to express the average (Tl n (r)} in terms of (3 n (r,t) is the high- 
temperature approximation j3F n {t) <C 1 or Hu) n -C kT. As a rule, this approxi- 
mation is well fulfilled in the NMR experiment. Making use of high-temperature 
expansion in Ea. (|72|) and taking into account that in this approximation 



exp(-S(t,0)) « — ( 1 - / d 6 rp n (f,t)H n (f f ) ) p L 
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we get 



'' '" j ' ' - £ 7T- D ^^7f-Mf, t) - (Mr) - P)R(f) (76) 



dt ' <9a; u 



or in a different form 

df3 n (r) 



Of 

Here D(f) is the diffusion coefficient 



D(r)Ap n (r) - (/3„(f) - /3)i?(r). (77) 



2h 2 u J lN(r)J- O0 J Tn(F-f 



N(r) = ^ fc 5(r — r*fc) being the nuclear spin density. The quantity i?(r) > has 
the following form: 

R(f) = /° e-^ / d 3 ri ^(^(r)^(^t l))L 

V ; h 2 ujlN(r) J-co l J Tn(I*y V 7 

Here the symbol (■ ■ -)l = Tr(. . . p^) implies the average over the equilibrium 
ensemble for lattice degrees of freedom. 



4.2 Host Nuclear Spin Diffusion in Dilute Alloys 

Spin diffusion is the transport of Zeeman energy or magnetization via the dipole- 
dipole interactions and it proved important both theoretically [SJ M I2U 1451 14"6*| 
and experimentally [01031 in diamagnetic solids. We consider here another class 
of substances, the dilute alloys [221113 IH3] The spin dynamics and relaxation 
of bulk metal nuclei by relatively dilute local moments in dilute alloys (e.g. 
Cu — Mn) was studied quite extensively, both theoretically H5J [H3 EH lEU 
El OH EU 021 OS] and experimentally |M| 0M 0M 0M HO %M QUI EHl EM 
ESI EH ESI EM EM 11011 11001 HUB- The description of spin relaxation in 
dilute alloys has certain specific features as compared with the homogeneous 
systems. For brevity we confine ourselves to the consideration of the bulk metal 
nuclei relaxation in dilute alloy. Due to the dipole-dipole interaction between a 
nuclear spin and an impurity spin, the relaxation rate may become nonuniform. 
It is more rapid for the spins that are close to impurity and is much slower for 
the distant nuclear spins. As a result, a nonuniform distribution in the bulk 
nuclear spin subsystem will occur and to describe spin relaxation consistently, 
the nuclear spin diffusion should be taken into account. 
The Hamiltonian for nuclear and electronic interacting spin subsystems is 

H = H n + H e + Hm + H n e + Hmc + Hdip (80) 

Here index n denotes the host nuclear spins, M denotes spin of the magnetic 
impurities, and e denotes the electron subsystem. In this section, when we refer 
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to the host nuclear spin subsystem H n we put 

H n = £ Iftuj n + £ A l3 (Jf J? - ij+jr). (81) 

i ijtj 

The Hamiltonian of electron subsystem is 

H e = E £kaa\ a a ka (82) 

A: it 

and 

H M = J2 huJ MS z m (83) 

m 

is the Hamiltonian of the impurity spins in the external magnetic field. The 
Hamiltonian of the interaction jllQ) of nuclear spins and the spin density a(Ri) 
of the conduction electrons is 

87T 

^ne — Jne / Jne — ~~2 3 (^^) 

i "> Inle 

where 

q q 

Interaction of the impurity spins S m and the spin density of the itinerant carriers 
is given by the spin-fermion [1111 1112] (sp — d(/))modcl Hamiltonian 

HMe = Jsd E S m a(R m ) (85) 
The last part of the total Hamiltonian (|80|l 

is the Hamiltonian of the dipole-dipole and pseudo-dipole interaction of nuclear 
and impurity spins. This interaction was described in detail in Refs. 94, 9*5l lll3| . 
The pseudo-dipolar interaction does not originate in crystalline anisotropy but 
in the tensor character of the dipolar interaction |95j . Their expression for the 
pseudo-dipolar interaction is 

Kn = YfiTj - 'irl 3 2 {kr7 3 ){T 3 r-; 3 )]B l0 (87) 

ij 

The Van Vleck Hamiltonian for a system with two magnetic ingredients |95U113| 
includes the term 

2 2 

■i9eV 2 ,5 \(to o„-2/ 



"I" ^ ' ( 3 I" Bi m ){IiS m ^ r ij (Ii r im)(S m T'im)} 

im lra 

^ t l~ B mn ){S m S n — 3r mn (Sirmn)(Sjr m n)}- (88) 
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The -B's represent the pseudo-dipolar interaction 

% = |(^ + ^)(l-3cos 2 %) 

. The later consists of three components of which we use in Eq. (|86[) the following 
one as the most essential |1 131 |9"S*] 

)) (89) 

im 

It was shown in [23 that for the large distance between the nuclear spin and 
the electron spin Bi m has the form 

B im « B ™^ + *») (90) 
(2fc F r; m )' ;i 

Thus, in structure, the coefficient Bi m is similar to the production of the contact 
potential and the spatial part of the RKKY interaction |114) . As a rule, the 
pseudo-dipolar interaction is less than the contact interaction. The estimations 
give B ~ 1/3 J ne for 205 Tl. It will be even more valid for copper since its mass 
is much less than for Tl. 

Now the expression for the Hamiltonian Hdi P can be rewritten as 

H dip = lnlM h js-U'SS^O- - 3 cos2 *<m) - (91) 

im %Ta 

3 

- sin0 im cos9 im [exp(-i(p im )lf6S^ + exp(i^ im )IrSS^)} 

cos(2k F r irn + 4> B ) 

{l + B 84 } 

Here we have introduced the mean field {S* n ) and the fluctuating part of the 
impurity spin, namely SS^ = — (S^). By substituting this definition of S* n 
into l|84|) rewritten in terms of the variable 5S m we obtain 

Hue = --3^— V(/> P )<5(7?, ; - f p ) = (92) 



(^(^-^-(a^-r;)))/?], 

where 

p kk' ss' 
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Now it is possible to write down explicitly the shift of the Zccman frequency to n 
in (|65|l due to the mean-field renormalization fli as 

fl^j^Mh^—il + B ^3 }(S )- 

Jne ^K 6 ^ ~ ^P)) = E ^raiSL) ~ Jne " ^ (93) 

This shift of the Zeeman frequency (fij <§C w n ) is the most essential for the 
evaluation of the coefficient of spin diffusion [33! I37l I7B] . 

4.3 Spin Diffusion Coefficient in Dilute Alloys 

Here we evaluate concrete expressions for the spin diffusion coefficient <|78|) for 
the dilute alloys system which is described by the Hamiltonian (|80|l . Consider 
again the approximate equation l|76|) where the diffusion coefficient can be writ- 
ten as 

D*" « J2 ~ ~ r?) exp[-(O r - (cj d ) 2 } (94) 

n v 71 " i 

In the derivation of the above expression, to permit explicit calculations, the 
Gaussian approximation for the nuclear spin correlation function was used ( see 
Appendix B). From equations l|76|) and l|94|) it follows that in the process of the 
longitudinal nuclear spin relaxation, which is a function of position, there is a 
possibility to transport the nuclear magnetization ( i.e. excess of nuclear spin 
density) due to the dipole-dipole interaction. It is clearly seen that the nuclei 
themselves do not move in the spin diffusion process. There is diffusion of the 
excess of the projection of the nuclear spin only. 

To proceed further, consider the case when the concentration of the impurity 
spins is very low. In this case, for one impurity spin there is a big number of 
host nuclear spins which interact with it. In other words, this case corresponds 
to the effective single-impurity situation. Thus, we can place one impurity 
spin to the origin of the coordinate frame (0,0,0). The vector r in Eq. l|94|l is 
then counted from this position. For a simple cubic crystalline system with the 
inversion center the symmetric tensor D^ v {r) is reduced to the scalar D{r). The 
coefficient D{r) decreases with decreasing the distance r when r is small. This 
is related with the fact that Zeeman nuclear frequencies of the nuclei, which are 
close to the impurity, have substantially different values due to the influence 
of the local magnetic fields induced by the impurity spin. This circumstance 
hinders the flip-flop (Q — ujm — Wn) transitions of neighboring nuclei since this 
transition does not conserve the total Zeeman energy of nuclear spins. ( Let 
us remind that if we suppose that the spins S are completely polarized and 
the nuclear spins I are completely unpolarized, then the dipolar interaction 
permits simultaneous reversals of S and / in the opposite directions, or flip- 
flops, and also reversals in the same direction which is usually called flip-flips 
with fl — u>m + ^n). In expression Ij94(l this tendency is described by the 
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exponential factor. This exponential factor leads to the appearance of the so- 
called " diffusion barrier" around each impurity. Inside this diffusion barrier the 
diffusion of nuclear spin is hindered strongly j^Bl ES| ■ 

It can be seen that for the large distance from the impurity the frequency 
difference in Eq.JMJ) behaves as (fl r — fi;) <C u>d, where LOd « 6j 2 fra~ 3 is the 
dipolar line-width and D(r) does not depend on r. In the opposite case, of 
small distance scale ( near impurity ) the frequency difference is big and the 
coefficient D(r) decreases quickly with the distance to the impurity. Thus, it is 
convenient to introduce the effective radius of the diffusion barrier <5, namely, a 
distance from the impurity for which the following definition holds: 

[d ifr>6 
DW = {» tfP <, < 95 > 

The constant D is equal to D = u)d/3h -\Z^^2A kl r kl . 

Let us estimate the "size" of the diffusion barrier. Consider two neighboring 
nuclei which take up a position along the radius from the impurity. The distance 
between them is equal to the lattice constant a. In this case, the frequency shift 
is equal to (fig — Qs+a) ~ ^>d and S w a^/[jM lln(S z )] ( see Appendix B). 
Consider again the approximate equation (|77|) taking into account the diffusion 
barrier approximation 1)95(1. It can be rewritten in the form 



dt 

where 



D A/3 n (f ; t) - (j3 n (r, t) - 0) (R + J2 X (0 + R 2 (r) ) , (96) 



2 J 2 f°° 
^ = i2fEE«H^' du>f{u-u n )Gl k , pp ,{u), (97) 



kk' pp' 



G°kk'pp'(^)= dtexp(ituj)(al L a kn a^(t)a p 'i(t)), (98) 



Rl{ff) = -TT^ll / - Un)Re(r k Tpk>Gt k , m (")^), (99) 

71 kk' rn 



and 



G\ k , m {u) = / dte^{itw){a^a kli S z m {t)), (100) 



q i f°° 

i? 2 (f) = rrj- ]T / duJ f^ ~ " n )Gl, lm (u;)Y m , (101) 

2(7«7M?i) 2 2tt ^ 

poo 

G* ror » = / dtexp(itu)(SS* m 6S* m (t)), (102) 



_ | Bcos(2k F \r-r m \ + B ) n 2 sin 2 6> rm cos 2 fl rm 
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Here the function f(uj — uj n ) is the NMR line-shape. The line-shape of the NMR 
spectrum |115| arises from the variation of the local field at a given nucleus 
because of the interaction with nearby neighbors. The inhomogeneity of the 
applied magnetic field may also increase the width of the line. 
The contribution of the factor Rq 1 leads to the generalized Korringa relaxation 
rate Q32] 

^^[f^x^lMo)! 2 }] 2 (104) 

Korringa |116| calculated the spin-lattice relaxation time 7\ in metals and 
showed that T\ should be inversely proportional to temperature and should 
be related to the Knight shift ( see also Ref. |117j ). Korringa nuclear spin- 
lattice relaxation occurs in a metal through the nucleus-electron interaction of 
contact type |116j 

-j(\ le \ns)( ln ni)\^ A (o)\ 2 . (105) 

The quantity R\ is determined by the correlation of the electron and impurity 
spins and is highly anisotropic. 

The quantity R2 is related to the scattering of nuclear spins on the fluctuations 
of impurity spins. The last contribution is the most essential factor in the 
present context. This is related to the fact that the main characteristic features 
of the problem under consideration clearly manifest itself in the isotropic case 
which is considered in the majority of works. In the isotropic case R\ = and 
the contribution of R2 can be expressed as 



C = -2— / duf{u> - u n )G l mm (cu) (107) 

5(7«7A/ft) 27r J-00 

Nevertheless, even after simplifications described above, a solution of the diffu- 
sion equation is still a complicated problem. The main difficulty is the presence 
of the highly oscillating factor cos(2/cp|r — r m \ + 4>b)- The role of this oscillat- 
ing factor can be taken into account entirely by numerical calculations. For a 
qualitative rough estimation we consider the simplified case when B ss 0. Then 
we can proceed following the method of calculation of Ref. 76 . According to 
these calculations (7S| we find 

— = (Ro)' 1 + AttDNF (108) 

Here N is the number of impurities and the quantity F has the form 



0.76 if b > 5 

l/3(b/5) 3 b if b<5 



p ^ ^ ( io9 ) 



where b = ty(C/D). 

It is clear from Eqs. l|108jl and H109JI that the behavior of the relaxation time 
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and its value depend strongly on the interrelation of b which is determined by 
the correlation function G l mm (uj) and of S which is determined by (S z ), as well 
as on the temperature for each concrete alloy. Thus, the problem of description 
of spin-lattice relaxation in dilute metallic alloys was reduced to the problem of 
calculation of the value of F. When 5 <C b the diffusion barrier is nonessential. 
In the opposite case, when b < S, the diffusion barrier is essential and leads 
to the slowing down of the relaxation process. In other words, the distance b 
determines the scale up to which the nuclear spin relaxation is effective. Finally, 
let us note that the order of value of time which is necessary to transmit the 
magnetic moment to the distance r in a solid is equal to td — r 2 /D; for r = 10~ 6 
cm it gives the value td — 1 sec. 

5 CONCLUDING REMARKS 

In the present paper, we have given a complementary method for obtaining 
the rate and relaxation equations of nuclear spin system in solids. The main 
tool in this approach is the use of the method of nonequilibrium statistical 
operator [HO] - We have presented a theory of spin relaxation which allows us 
to derive general equations of spin dynamics. In addition, our theory allows 
us to take into account the effects of spin diffusion in a very straightforward 
manner. The calculations were kept general by restricting the form of spin- 
lattice Hamiltonian as little as possible. It has permitted us to perform the 
derivation under more general conditions and explicitly demonstrate some key 
features of irreversible processes in solids. 

It was shown that the spin systems provide a useful proving ground for applying 
the sophisticated methods of statistical thermodynamics. The method used is 
capable of systematic improvement and gives a deeper insight into the meaning 
of the spin relaxation processes in solids. We have shown that the transport 
of nuclear spin energy in a lattice of paramagnetic spins with magnetic dipolar 
interaction plays an important role in relaxation processes in solids. To test the 
general formalism presented here, an example of a dilute metallic alloy system 
was considered to demonstrate the usefulness of the equations derived. 
In summary, the present paper examines the relaxation dynamics of a spin 
system. It continues the investigation presented in the previous work into the 
use of statistical mechanical methods for systems that are in contact with a 
thermal bath. We used the method of the nonequilibrium statistical operator 
developed by D. N. Zubarev. In the present paper, we have developed the 
application of this method to the spin-relaxation problem, so that some useful 
results may be obtained from it. The calculation presented in this paper can 
be said to show that the NSO method has provided a compact and efficient 
tool for description of the spin relaxation dynamics. In this respect, the present 
treatment may be regarded as a complement of the Buishvili and Zubarev |28| 
seminal treatment. 

Though the analysis of this paper concentrates on the nuclear spin systems in 
solids, the extension to other spin systems, e.g., paramagnetic electron spin 
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system, is straightforward. The other important task is to examine the effects 
of a periodically time-dependent field on the long-time behavior of an otherwise 
isolated system of many coupled spins. This question is a part of a more general 
problem of the evolution of a complex system in an external field, especially in 
an intense external field. We hope that the methods here developed may be 
applied in these cases with the suitable modifications. 

A EVOLUTION OF A SYSTEM IN AN AL- 
TERNATING EXTERNAL FIELD 

In section (|2.2() we wrote the kinetic and evolution equations in the approach of 
the nonequilibrium statistical operator. In this Appendix we show briefly the 
derivation of the same equations in the presence of alternating external field. 
This problem is essential for the nuclear and electron spin resonance. Both 
nuclear and electron spins have associated magnetic dipole moments which can 
absorb radiation, usually at radio or microwave frequencies. 
We consider the many-particle system with the Hamiltonian 

H = H 1 +H 2 + V + H f {t), (110) 

where 

H 1 = Y,E a a) a a a (111) 

a 

is the single-particle second-quantized Hamiltonian of the quasiparticles with 
energies E a . This term corresponds to the kinetic energy of noninteracting 
particles 

N P 2 N h 2 

i=l i=l 

The index a = (fc, s) denotes the momentum and spin 

fa( x ) = Vfc(r)A(s — cr ) = exp(i/cr)A(s — a)/y/v 

E a = (alHila) 

1 f ft 2 h 2 h 2 

(k\H(l)\k') = - / d 3 rexp(ikr)( V 2 )exp(ifc'r) = A(k-k') 

v J 2m 2m 

and 

^ = 5> Q/3 4a/3, *o!9 = *J Q . (H2) 

Operator V is the operator of the interaction between the small subsystem and 
the thermal bath, and H2 is the Hamiltonian of the thermal bath which we do 
not write explicitly. The quantities $ Q/ 3 are the operators acting on the thermal 
bath variables with the properties ($ a pY = ®% a ; ®% a = ® a p- 
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The interaction of the system with the external time dependent alternating field 
is described by the operator 



#/(*)= X>«/j(*)4«0> (H3) 

For purposes of calculation, it is convenient to rewrite Hamiltonian Hf{t) in a 
somewhat different form 

Hf(t) = -J2 T (^MaUp, (H4) 



where 



and 



v 

a, (3 



h a p{t) = -T{a,0,t) 
v 



T = Y J T(r l ,t): T(p) = J d 3 r exp(ipr)T{r,t) 

1—1 

(k\T(r,t)\k') = - [ d 3 reMi(k-k')r)T{r 7 t) = -T(k-k',t) 



We are interested in the kinetic stage of the nonequilibrium process in the sys- 
tem weakly coupled to the thermal bath. Therefore, we assume that the state of 
this system is determined completely by the set of averages {P a p) — (a^ap) and 
the state of the thermal bath by (H2), where (. . .) denotes the statistical average 
with the nonequilibrium statistical operator, which will be defined below. 
Following Pokrowsky's calculations we can write down the nonequilibrium sta- 
tistical operator in the following form: 

p(t)=Q- 1 cxp(-L(t)), (115) 

where 

a/3 

(116) 

The notation TL2 denotes TL2 = H2 — P2N2 where 112 is the chemical potential of 
the medium (thermal bath) and J2 = Ti.2(ti). In this equation, the time depen- 
dence of the operators in the right-hand side differs from the time dependence 
in Ea. (|25() . Consider this question in detail. The Heisenberg representation 

—iHt 

H{t) = W(t)HU(t); U(t) = exp(^^) (117) 
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in the presence of the external field H = Hq + H ext takes the form 

A{tv,t) = U^{t + h;t)AU{t + ti;t) (118) 

where 

i r t+tl 

U(t + t i; t)=Tcxp(-- H ext {T)dr) (119) 

" Jt 

We now generalize the evolution equations to the case in which the external 
field is present. We have 

P a /3 = -^\ P ^H] = 
^{E p - E a )P afi + VW - h va {t)P vl3 ) + ^[P a p, V] (120) 

Then we can write down the balance equation 

Ji + -h = If, (121) 

where 

Ji=H 1 = ^[Hi,H] = ^{[HuV] + [H U H^]) (122) 

and 

// = Ji + J 2 = ^ ]T{(£ Q - E (j )h afj {t) + h af} (t)[P af ,, V}} (123) 

a/3 

The last term describes the work of the external field. 

The parameters F a p(i) are determined from the condition {P a f}) = (Pap)q- The 
quasi-cquilibrium statistical operator p q has the form 

P q =Pip2, (124) 

where 

pi = Q^ 1 exp - P a pF afj (t) ; 

V a/3 J 

Qi = Tr cxp \^-J2 

P a f)F a p{t)^ (125) 

p 2 = Q2 1 cxp {-P(H 2 - P2N2)) ; 

Q 2 =Trexp(-/?(ff 2 - M2 7V 2 )) (126) 

Thus, we can write 

-^--<[P a/3 ,if]>- 

(127) 

l(^ - B Q )<p a/J ) + ^ E (( p -)vw - ^w(^)) + v]) 
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To calculate explicitly the r.h.s. of Ea. (|127l) we use the formula 

exp(-A - B) = exp(-A) - / exp(-A)(exp(-AT)Bexp(AT)d,T) (128) 

Jo 

p~{l- {exp(-AT)Bexp{Ar) ~ (cxp(- At) Bcxp{AT)) A )dT)}p{A), (129) 
Jo 

where 

exp(-A) 



P{A) 



Trexp(-A) 



a/3 



r° i 

00 fit; 

Making the same expansion procedure as described in section l?~2l wc find 

(130) 

+ 4 Y, i{P av )hpv{t) - K a (t)(P v(3 )) + 4 / d\{[P afi , V]e- XA Ve XA ) q + 
in ^ in Jo 

(ill) Jo J-oo 

where 

X a , p ,(t) = F a ,fj,{t) - @[8 a ,p,E al + h a ,p,(t)} (131) 
It can be rewritten as 

d{Pap) ~Ue p -E a )(P aP )+ (132) 



dt ih 

v a' f}' 

where the generalized relaxation matrix is given by 

_ P f 1 



f dX([P aP ,V}e~ XA Ve XA ) q + 
Jo 

(133) 



1 '° 



d\ / dt 1 e Etl <[P a ^/,^]e- AA [P a ^(ti ! t),y(t 1 )]e A ^) 9 X Q , /3 ,(t + t 1 ) 



DO 



W Jo 

Equation l|132|) gives the generalization of the rate equation (|34|) of the Redfield- 
type for the case of the external alternating field. A more detailed investigation 
of this equation and the problem of the evolution of a system in an external 
field will be carried out separately. 
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B CORRELATION FUNCTIONS AND GAUS- 
SIAN APPROXIMATION 

The eigenvalues of the Hamiltonian (|80|l correspond to well-defined values of 
y~]. If = I z = to. Their energy is the sum of a Zeeman energy mhuj n and 
a spin-spin energy. A stochastic-theoretical treatment of the spin relaxation 
phenomena is a useful complementary approach to the consideration of spin 
evolution [1181 1119| . By a stochastic theory it is termed ordinary that kind of 
theoretical treatment of the problem in which one assumes the random nature 
of the forces acting on a system. The phenomenon of spin relaxation can be 
properly interpreted as some stochastic process of spin motion. This stochastic 
process is determined by the equation of motion of the spin variable. It was 
formulated |331 1118| plausibly that a Gaussian random process may be well 
applied for the evolution of the magnetization in the presence of a static external 
field 

—jl=jflx(h + h), (134) 

where 7 denotes the gyromagnetic ratio, ho a static external field, and h the fluc- 
tuating internal field due to the magnetic moments in the surrounding medium. 
The effect of the fluctuating internal field h is to cause nuclear spin transitions 
governed by the selection rule Am = ±1. If the Zeeman splitting is small, 
i.e., fiuj n <C kT, then the transition probability for a Am = ±1 transition will 
be proportional to the Fourier transform of correlation functions of the form 
(h+(t)h-(t')), (h-(t)h+(t')), (h z (t)h z (t')). If we assume the process of h(t) to 
be a Gaussian random process, the problem becomes more easy tractable. From 
this viewpoint it is reasonable to assume that the equation of spin motion in- 
volves the local fluctuating magnetic field whose process is assumed to be a 
Gaussian random process |118| . 

The Gaussian or normal probability distribution law is the limit of the binomial 
distribution 

P(m) = C™p m (l-p) n 

in the limit of large n andpn (n — > 00). Here n is the repetition of an experiment, 
p is the probability of success, and C™ = n\/m\(n—m)\. The normal probability 
distribution has the form 



P(m) = 




(135) 



where a = ^Jnp(l — p) is a measure of the width of the distribution. It is 
clear that the Gaussian distribution results when an experiment with a finite 
probability of success is repeated a very large number of times. The Gaussian 
random process is a random process (with discrete or continuous time) which 
has the normal (Gaussian) probability distribution law for any group of values 
of the process. The Gaussian random process is determined completely by its 
average value and correlation function. Thus, the description of the class of 
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Gaussian processes is reduced to the determination of the possible form of the 
corresponding correlation functions. 

Consider now an isotropic distribution of nuclei and rewrite the production of 
the current operator in Ea. (|78|l in explicit form 

J(r)J(fl,ti) = ^pX! E A kiA mn r k ir mn (136) 

dfS(r- fk)8(fi - )[Tr(Z 2 ) 2 ]- 1 Tr7+/ ; -7+(t 1 )7-(t 1 ) 

To proceed further, the form of the correlation function of nuclear spins in the 
above expression must be determined. In the theory of NMR the reasonable 
assumption is that this correlation function can be represented in an intuitively 
understandable way as [231 1118) 

TrI+I-I+{t)I-(t) cx Tr(I + I-) 2 f(t)S kn S lm eMj(^i - «fc)] = 

-J(t)5 kn 5 lm exp[~(fi, - fife)] (137) 
4 ft 

Then the diffusion coefficient D{r) (|78|l takes the form 

D ^ = ^2T77T E A tl r kAr - fk) / e*dtf(t) exp[-(Sl k - fi,)] 

= A rl r rl / dtf(t)eX P [-(n r n,)]. (138) 

8 ft , J-oo ft 



The method of moments gives that f(t) is close to the normal probability dis- 
tribution 

m = AeM-^); hW d = ^-^ (139) 
The constant A can be determined from the condition 



/OO / r) 

dtf(t) = l = A.—; A 
-oo V w d 



(140) 



Thus, we obtain 



/(t)= e«p(-^) (141) 



/2tt 

For the diffusion coefficient Ijl38(l we find 

D(f) « ^^^exp[-(n r - ^) 2 /4 M 2 ]. (142) 
ft V 71 " j 

In the case when r is close to Z the frequency difference (f2 r — f2;) 3> and 
Z3(r) — > 0. In the opposite case, when (f2 r — f2;) <C o-id the diffusion coefficient is 
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nearly constant D(r) ~ D. Thus, we traced back to the notion of the diffusion 
barrier S (|95() . Consider two neighboring nuclei along the radius from the im- 
purity. The distance between them is equal to the lattice constant a. For this 
case the frequency shift is equal to 

where u>d ~ Qln^ a ~ 3 - 

Consider this constraint more carefully. We have 

t icrz\ r 1 + Bcos(Sk F + 4>b) 1 + Bcos((5 + a)k F + </> B ) 

wuHSH ^ } = 

M „„ ,1 + Bcos(5k F + 6 B ) 1 — Bsin(8k F + cj) B )ak F 
1u1mHS z ){ p ^ + 

B cos(5k F + <Pb) r l + Bcos(5k F + (j) B ) ^ 
+ p + 3 ^ a} = 

-jplnlMHS ){3 

B sin(<Sfc F + <j} B )k F } = 6jlha- 3 (143) 
For the rough estimation we omit the cos and sin contributions. Then we obtain 



6 7 >- 3 = 7 „ 7 mMS z }^; S = a^[^(S')] (144) 
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